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Abstract 



For the first time, the nuclear attenuation of three hadron systems 
is studied in neutrino-induced reactions using the data obtained with 
SKAT bubble chamber. The strongest attenuation (R3 ~ 0.6) is ob- 
served for a system carrying an overwhelming fraction of the current 
quark energy, as well as for a system with the smallest effective mass. 
An indication is obtained that the correlation effects in the nuclear at- 
tenuation play only a minor role. The experimental data are compared 
with predictions of the quark string fragmentation model. 



1. INTRODUCTION 



The investigation of the multihadron leptoproduction on nuclear targets 
can provide a valuable information on the space-time structure of the quark 
string fragmentation and the hadron formation ([1-4] and references therein). 
At present the experimental data on the double hadron production are avail- 
able in electronuclear [5] and neutrinonuclear [6] interactions. This work is 
devoted to the study, for the first time, of the nuclear effects in the neu- 
trinoproduction of more complex, three hadron systems. In Section 2 the 
experimental procedure is described. The experimental data are presented 
and compared with theoretical predictions in Section 3. The results are sum- 
marized in Section 4. 

2. EXPERIMENTAL PROCEDURE 

The experiment was performed with SKAT bubble [7] chamber exposed to 
a wideband neutrino beam obtained with a 70 GeV primary protons from 
the Serpukhov accelerator. The chamber was filled with a propane-freon 
mixture containing 87 vol% propane (CaHg) and 13 vol% freon (CFsBr) with 
the percentage of nuclei H:C:F:Br = 67.9:26.8:4.0:1.3%. A 20 kG uniform 
magnetic field was provided within the operating chamber volume. 

Charged current interactions containing a negative muon with momentum 
> 0.5 GeV/c were selected. Other negatively charged particles were con- 
sidered to be n~ mesons. Protons with momentum below 0.6 GeV/c and a 
fraction of protons with momentum 0.6-0.85 GeV/c were identified by their 
stopping in the chamber. To non-identified positively charged particles the 
pion mass was assigned. Events with errors in measuring the momenta of all 
charged secondaries and photons less than 60% and 100%, respectively, were 
selected. Each event was given a weight to correct for the fraction of events 
excluded due to improperly reconstruction. More details concerning the ex- 
perimental procedure, in particular, the estimation of the neutrino energy 



E v and the reconstruction of ir° — > 27 decays, can be found in our previous 
publications [8,9,10]. In this work, 77 combinations with the effective mass 
m 77 closest to 135 MeV/c 2 (and enclosed in the range of 105 < m 77 < 165 
MeV/c 2 ) were taken, requiring the each detected 7 to be included in no more 
than one 77 combination. 

The events with 3 < E v < 30 GeV were accepted provided that the recon- 
structed mass W of the hadronic system exceeds 2 GeV and y = v j ' E v < 0.95, 
v being the energy transferred to the hadronic system. No restriction was 
imposed on the transfer momentum squared Q 2 . The number of accepted 
events was 4499 (5683 weighted events). The contamination from the neutral 
current (NC) interactions to the selected event sample was estimated to be 
(2.6±0.2)%. The mean values of the kinematic variables were (E u ) = 10.0 
GeV, (W) = 2.9 GeV, (W 2 ) = 9.1 GeV 2 , (Q 2 ) = 2.7 (GeV/c) 2 and (u) = 
5.8 GeV. 

The nuclear attenuation effects studied below are inferred with the help of 
a comparison of the trihadron characteristics in two event subsamples: the 
nuclear subsample (Ba) and quasinucleon subsample (Bjf) composed using 
a number of topological and kinematical criteria [SI |9]. The Bn subsample 
includes events with no indication of the nuclear disintegration or a secondary 
intranuclear interaction: the total charge of secondary hadrons is required to 
be q = +1 (for the subsample B n of interactions with a neutron) or q = +2 
(for the subsample B p of interactions with a proton), while the number of 
recorded baryons (these included identified protons and A -hyperons, along 
with neutrons that suffered a secondary interaction in the chamber) was 
forbidden to exceed unity, baryons flying in the backward hemisphere being 
required to be absent among them. Moreover, a constraint was imposed 
on the effective target mass M t < 1.2 GeV/c 2 , the M t being defined as 
M t = Tj(Ei —pf) where the summation is performed over the energies Ei and 
the longitudinal momenta pf (along the neutrino direction) of all recorded 
secondary particles. Events that did not satisfy aforementioned criteria were 



included in the subsample B$ of cascade events. As a result, the weighted 
numbers of the B p , B n and B$ proved to be 1249, 1304 and 3129. 

The validity of the selection of the B p and B n subsamples was verified 
[SJ dU 12] by comparison of a large number of the multiplicity and spectral 
characteristics of hadrons in these subsamples with the available data ob- 
tained on hydrogen and deuterium targets, resulting in a satisfactory agree- 
ment and giving a sufficient ground to conclude that the B p and B n sub- 
samples may contain only a minor contamination from events where the 
secondary intranuclear interactions deteriorate the characteristics of the pri- 
mary hadrons. We checked the sensitivity of the extracted data (presented in 
the next section) relative to the choice of the boundary value of M t , varying 
the latter from 1.1 to 1.4 GeV/c 2 , e.g. changing the event numbers in the 
Bp, B n and B$ subsamples and continuously increasing the contamination 
from secondary interactions in the subsamples B p and B n (see Fig. 1 in [TTj). 
We found that the variation of the dihadron characteristics have been on an 
average fivefold smaller as compared to the statistical errors. 

As it follows from the composition of the propane-freon mixture (see above) 
about 36% of sub-sample B p is contributed by interactions with free hydrogen 
(at the ratio of vn and vp CC cross sections being equal to 2). Weighting 
the quasiproton events with a factor of 0.64, one can compose a pure nuclear 
subsample Ba = B s + B n + 0.64£> p and a quasinucleon subsample B N = 
B n + 0.64i?p. As it follows from the percentage of C, F and Br nuclei in 
the propan-freon mixture the mean atomic weight of the composite target 
is equal to A = 27.6 for neutrinonuclear interactions. However, since the 
nuclear attenuation effects are under consideration, it seems more relevant 
to introduce an effective nucleus A e s in which the probability of the hadron 
absorption is the same as that obtained from averaging over the target nuclei. 
The calculations (see [13] for details) result in 19< A cS <23, depending on 
the momentum of the neutrinoproduced hadron. 



3. THE RESULTS 



Below we will consider hadrons produced in the forward hemisphere in the 
hadronic c.m.s. (i.e. in the region of xf > 0, xp being the Feynman variable), 
because in this region the nuclear attenuation effects for hadrons dominate 
over multiplication effects caused by secondary intranuclear inelastic inter- 
actions, observable in the region of xf < and, to a much smaller extent, 
even in the region of < x F < 0.1 [9,11]. Below we will imply two different 
cuts on xp- xp > and xp > 0.1. 

For a given kinematic variable v, characterizing a trihadron, the nuclear 
attenuation is measured by the ratio Rz(v) = nf(v)/n^ (v), where n^{v) and 
n% (v) are the trihadron differential mean multiplicities (yields) in the nuclear 
and quasinucleon subsamples, respectively. 

Figure 1 (top panel) shows the dependence of R% on the maximal angle 
i? max between two hadrons entering into the three-hadron system. The data 
exhibit a strengthening attenuation with decreasing $ max which can be caused 
by secondary intranuclear interactions, mainly elastic for the case of xp > 0.1 
and both elastic and inelastic for the case of xp > 0. The elastic scattering 
of hadrons also leads to the increasing effective mass m e fr of the trihadron 
and hence to a stronger attenuation at smaller m c fj than at larger m e fr, as it 
is seen from the middle panel of Figure 1 for x F > 0.1. On the other hand, 
the inelastic interactions (playing more significant role at the cut x F > 0) 
lead to energy losses of hadrons which partly compensate the effective mass 
increasing caused by the widening of the angles between scattered hadrons, 
as it can be concluded from the middle panel of Figure 1 where a rather flat 
m e ff - dependence is observed for the case of xp > 0. 

The bottom panel of Figure 1 shows the dependence of the ratio R^{pt) 
on the trihadron transverse momentum p t (defined relative to the current 
quark direction). A trend of increasing R^{pt) with increasing p t is observed, 
being a consequence of secondary interactions of neutrinoproduced hadrons. 



It should be noted, that a characteristic feature of the data presented in 
Figure 1 is that, in general, the ratio Rs(v) at xp > exceeds that for the 
case of xp > 0.1, indicating, as it was already mentioned above, that a non- 
negligible fraction of low xp hadrons (with < x F < 0.1) originates from 
inelastic intranuclear interactions of higher - xp hadrons. 

Figure 2 presents the dependence of the ratio i? 3 on the collective variable 
Z = Yl z k, where Zk = E^jv with being the energy of the k- th hadron 
and the sum is over three hadrons. For comparison, we also plotted the 
ratio R 2 {z) versus Z = z\ + z 2 for two-hadron systems and the ratio Ri(z) 
versus the z value for single hadrons. As it is seen, both functions R±, R 2 
and i?3 continuously decrease with the increasing argument. It worths to 
note, that the fastest single hadron (with z > 0.85) absorbs stronger than 
the fastest hadron systems (with the same summary Z) composed of two or 
three hadrons of smaller z's. 

It is also interesting to compare the ratio R\{z) for single hadron with the 
ratios R 2 (z) and R${z) where the argument refers to a hadron (the "trigger" 
hadron) accompanied by, respectively, one and two hadrons of arbitrary val- 
ues of the variable z. A decreasing n- dependence of R n (z), seen in Figure 
3 and, more clearly, in Figure 4 and 5, is related to the additional attenu- 
ation of the accompanying hadron(s). The n- dependencies of R n (z) plot- 
ted in Figure 4 (at Xp > 0) and 5 (at x-p > 0.1) can be approximated as 
R n (z) = pi(z)(l — ifi)™" 1 , where pi(z) is the attenuation ratio for the 'trig- 
ger' hadron, while w\ is the mean absorption probability of a single hadron 
averaged over hadrons accompanying the 'trigger' one. As it is seen from 
the fitted values of pi(z) and w\ quoted in Figures 4 and 5, the ratio pi(z) 
is maximal for the slowest trigger hadron (even exceeding unity at xp > 
due to the nuclear enhancement effects). At higher z- values, the variation 
of pi(z) is not significant, while the values of w 1 are around w 1 ~ 0.05 at 
x F > and around W\ ~ 0.12 at xp > 0.1. These W\ values are averaged 
over a wide range of the z variable of accompanying hadrons. A more definite 



estimation for the absorption probability can be inferred from the n- depen- 
dence of Rn(z) at a narrow z- range of all involved hadrons, for example, 
at 0.1 < z < 0.33 (see Figure 6). Under assumption that hadrons absorb 
independently with a mean probability w^ s , the n- dependence of R n (z) can 
be expressed as R n (z) = (1 — w a b s ) n - The data at x F > are described by 
this dependence quit well (with the fit value w a b s = 0.060 ± 0.006), being 
consistent with the assumption that the correlated absorption effects do not 
play a prominent role (Figure 6, the left panel). The description of the data 
at xf > 0.1 (with the fit value w abs = 0.077 ± 0.014) is somewhat worse, due 
to a possible (statistically not provided) deviation of the R 3 value from the 
general trend of the n- dependence of R n . It can be a faint manifestation 
of the mutual screening effect (suggested in [4]) in the system of hadrons 
produced at close impact parameters. Statistically more provided data are 
needed to check this assumption. 

It should be also noted, that the observed difference between w± and w a b s 
(quoted in Figures 5 and 6 for the case of xp > 0.1) can be caused by the 
fact that the latter concerns the hadrons with 0.1 < z < 0.33 at which the 
hadron formation length is expected to be maximal (hence reducing their 
absorption), while the former concerns the hadrons which can carry larger 
z's at which the formation length is comparatively shorter [14]. 

We undertook an attempt to describe our data in the framework of the 
Two-Scale Model (TSM) [15] (see also [4] and references therein) which was 
already applied for description of our data on the nuclear attenuation of neu- 
trinoproduced two-hadron systems [6]. The details of the TSM application 
to three-hadron systems can be found in [16]. The TSM contains four free 
parameters: the quark string tension k and the following string-nucleon inter- 
action cross sections, namely, o q for the initial string stretched between the 
struck quark and the target nucleon remnant, o s for the intermediate string 
stretched between the struck quark and a created antiquark (which becomes 
a valence one for the hadron being looked at) and a h for the formed colorless 



system with quantum number and valence content of the final hadron. 

As it was shown in [6], a reasonable agreement with the data on the two- 
hadron neutrinoproduction can be reached at the following set of the model 
parameters: k = 0.8 GeV/fm. a q = 0, a s = 10 mb and a h = 20 mb. This 
set was used to get the TSM predictions presented below. The predictions 
concern pions, pion pairs and pion triplets. In present calculations, only the 
z variable of produced pions is considered, while the xp variable does not 
appear explicitly in the model. Therefore, the same predictions of the model 
are compared below, in Figs. 2, 3, 4, 5 and 6, with experimental data inferred 
both at x F > and x F > 0.1. 

As it is seen, the description of the data at xp > is noticeably worse, 
probably, due to the aforementioned effects of intranuclear interactions (not 
inserted into the model) leading to a significant enhancement of the yield of 
\ow-z hadrons (especially in the region of z < 0.2) and, as a consequence, 
the yields of dihadrons and trihadrons with comparatively small values of 
Z < 0.3 -h 0.4. At a more severe restriction on xf > 0.1, the said discrepancy 
is much smaller, and the model provides a reasonable description of the 
data at Z > 0.2, as it can be seen from Figs. 2 and 3. An exception 
is R\(z) for the fastest single hadron with z > 0.85 for which the model 
badly underestimates the nuclear attenuation. On the other hand, the model 
predictions are compatible with the data at the the same summary Z-range 
[Z ~ 0.9) for hadron systems composed of two or three hadrons with smaller 
individual z's (see Fig. 2). The TSM describes quite well the n-dependence 
of the ratio R n (z) at ^-values of the trigger hadron above z = 0.22, as it 
can be seen from Fig. 3 (the right panel) and Fig. 5 (the middle and right 
panels). It should be also noted, that the TSM predictions for the slope in 
the n-dependence of R^z) (the dashed lines in Fig. 5) at z > 0.22 are rather 
close to those extracted under assumption of an independent absorption of 
neutrinoproduced hadrons (the solid lines in Fig. 5). The same can be said 
about the slopes in the n-dependence of R n for the case when all involved 



hadrons are enclosed in a narrow z-range of 0.1 < z < 0.33, as it is seen from 
Fig. 6 where the model predictions are plotted by dashed lines. The predicted 
values of R n , however, underestimate the experimental data, probably, due 
to the unsufficiently large values of the variable z (see the discussion above). 

4. SUMMARY 

The nuclear attenuation of three-charged hadron systems is studied for the 
first time in neutrino-induced reactions. 

The dependence of the attenuation factor R 3 on various kinematic variables 
of the trihadron is investigated. The strongest attenuation, with i? 3 ~ 0.6, 
is observed for a system of hadrons with x-p > 0.1 having the smallest ef- 
fective mass, as well as at largest values of the collective variable Z. The 
observed dependence of R n on the number n of involved hadrons enclosed 
in the narrow range of 0.1 < z < 0.33 is compatible with an assumption 
that the intranuclear interaction of a hadron occurs with a mean probability 
Wabs = 0.07 ± 0.02 and almost independently of accompanying hadrons. 

It is shown, that the Two-Scale Model of the quark string fragmentation 
provides a reasonable description of the nuclear attenuation for systems com- 
posed of relatively energetic (x F > 0.1 and z > 0.2) hadrons. 
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Figure 1: The dependence of R 3 on $ max (top panels), m cS (middle panels) 
and p t (bottom panels) for hadrons with x F > (left panels) and x F > 0.1 
(right panels). 
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Figure 2: The dependence of R n (n = 1, 2, 3) on Z for hadrons with a; F > 
(left panel) and x F > 0.1 (right panel) The curves are the TSM predictions 
(see text). 
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Figure 3: The dependence of R n (z) on the z variable of the "trigger" hadron 
at xp > (the left panel) and xp > 0.1 (the right panel) The curves are the 
TSM predictions (see text). 
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Figure 4: The n- dependence of R n {z) for different z- ranges of the "trigger" 
hadron at x-p > 0. The solid lines are the fit results, while the dashed lines 
are the TSM predictions (see the text). 
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Figure 5: The same as Fig. 4, but at x F > 0.1. 
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Figure 6: The dependence of R n (z) on the number n of involved hadrons 
with 0.1 < z < 0.33 and xp > (the left panel) and x-p > 0.1 (the right 
panel). The solid lines are the fit results, while the dashed lines are the TSM 
predictions (see the text). 



